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1 Introduction 



1.1 Motivation: A model problem 

Many physical problems can be modeled by partial differential equations. Let 
us consider for example the case of an elastic membrane ft, with fixed bound- 
ary T, subject to pressure forces /. The vertical membrane displacement is 
represented by a real valued function u, which solves the equation 

— Au(x) = f(x), x = (xi,x 2 ) G ft, (1.1) 



1 



where the Laplace operator A is defined, in two dimensions, by 



Au 



d 2 u d 2 



u 



dx 2 



dx\ ' 



As the membrane is glued to the curve T, u satisfies the condition 



u{x) = 0, x G r. 



(1.2) 



The system f ll.ll) - fll.2p is the homogeneous Dirichlet problem for the Laplace 
operator. It enters the more general framework of (linear) elliptic boundary 
value problems, which consist of a (linear) partial differential equation (in 
the example above, of order two: the highest order in the derivatives) inside 
an open set Q of the whole space M. N , satisfying some "elliptic" property, 
completed by (linear) conditions on the boundary T of i?, called "boundary 
conditions" . In the sequel, we only consider the linear case. 

Our aim is to answer the following questions: does this problem admit 
a solution? in which space? is this solution unique? does it depend con- 
tinuously on the given data /? In case of positive answers, we say that the 
problem is "well posed" in the Hadamard sense. But other questions can also 
be raised, such as the sign of the solution, for example, or its regularity. We 
give a full survey of linear elliptic problems in a bounded or in an exterior 
domain with a sufficiently smooth boundary and in the whole space. In the 
general theory of the elliptic problem we consider only smooth coefficients. 
We survey the standard theory, which can be found in the several well-known 
monographs of 1960s. The new trends in investigation of the elliptic prob- 
lems is to consider more general domains with nonsmooth boundaries and 
nonsmooth coefficients. On the other hand the regularity results for ellip- 
tic systems have not been improved during last thirty years. New trends 
also require employment of more general function spaces and more general 
functional background. 

The number of references (Further reading) is strictly limited here; there 
are only some of the most important publications. The basic facts can usually 
be found in more places and sometimes we do not mention the particular 
reference. To the very basic references throughout belong jF], [GT] . [DLJ, 
|Hoj . |LU] , |LM] , [RRj . [W| ; of course there are many others. 
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1.2 The method 



To answer the above questions, we generally use, for such elliptic problems, 
an approach based on what is called a "variational formulation" (see section 
3 below): the boundary value problem is first transformed into a variational 
problem of lower order, which is solved in a Hilbertian frame with help of 
the Lax-Milgram Theorem (based on the representation theorem). All ques- 
tions are then solved (existence, uniqueness, continuity in terms of the data, 
regularity). But this variational formalism does not necessary allow to treat 
all the situations and it is limited to the Hilbertian case. Other strategies 
can then be developed, based on a priori estimates and duality arguments 
for the existence problem, or maximum principle for the question of unic- 
ity Without forgetting the particular cases where an explicit Green kernel 
is computable (the Laplacean operator in the whole space case for example). 

Moreover, the study of linear elliptic equations is directly linked to the 
function spaces background. It is the reason why we first deal with Sobolev 
spaces — both of the integer and fractional order and we survey their basic 
properties, imbedding and trace theorems. We pay attention to the Riesz 
and Bessel potentials and we define weighted Sobolev spaces important in 
the context of unbounded opens. Second, we present the variational approach 
and the Lax-Milgram theorem as a key point to solve a large class of boundary 
value problems. We give examples: the Dirichlet and Neumann problems for 
the Poisson equation, the Newton problem for more general second order 
operators; we also investigate mixed boundary conditions and present an 
example of a problem of fourth order. Then, we briefly present the arguments 
for studying general elliptic problems and concentrate on second order elliptic 
problems; we recall the weak and strong maximum principle, formulate the 
Fredholm alternative and tackle the regularity questions. Moreover, we are 
interested in the existence and uniqueness of solution of the Laplace equation 
in the whole space and in exterior opens. Finally, we present some particular 
examples arising from physical problems, either in fluid mechanics (the Stokes 
system) or in elasticity. 

2 Sobolev and other type of spaces 

Throughout, Q C M. N will generally be an open subset of the iV-dimensional 
Euclidean space WL N . A domain will be an open and connected subset of 
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M. N . We shall use standard notations for the spaces L P {Q), C°°(Q) etc. and 
their norms. Let us agree that C k ' r (f2), k G N, r G (0, 1), denote the space 
of functions / in C fc (f2), whose derivatives D a f, a = (a±, . . . , ckjv) G of 
order |a| = X^ili = ^ are all r-H61der continuous. By Q in the notations 
for some of these spaces we mean that the functions have the corresponding 
property on Q and that they can be continuously extended to Q. 

Let us recall several fundamental concepts. The space T>(Q) of the test 
functions in Q consists of all infinitely differentiate tp with a compact sup- 
port in Q. A locally convex topology can be introduced here. The ele- 
ments of the dual space V{Q) are called the distributions. If / G L\ oc (Q) 
(i.e. / G L l (K) for all compact subset K of Q), then / is a regular dis- 
tribution; the duality is represented by j n f(x)tp(x) dx. If / G T>'(f2), we 
define the distributional or the weak derivative D a of / as the distribution 
(p I— > (— 1) '"'(/, D a (p). Plainly, if / G has "classical" partial derivatives 
in Ll oc , then it coincides with the corresponding weak derivative. 

If Q = M. N , it is sometimes more suitable to work with the tempered 
distributions. The role of T>(Q) is played by the space «S(R Ar ) of C°° functions 
with finite pseudonorms sup \D a f(x)\(l + \a\,k = 0,1,2,.... Recall 

that the Fourier transform JF maps S(M N ) into itself and the same is true 
for the space of the tempered distributions S'(M. N ). 

2.1 Sobolev spaces of positive order 

The Sobolev space W k ' p {Q), 1 < p < oo, k G N, is the space of all / G L p {0), 
whose weak derivatives up to order k are regular distributions belonging to 
L p (i7); in W k ' p (f2) we introduce the norm 




when p < oo and max| a |< fc sup ess xg ^ \D a f(x)\ if p = oo. The space W k,p (f2) 
is a Banach space, separable for p < oo, reflexive for 1 < p < oo; it is a 
Hilbert space for p — 2, more simply denoted H m (Q). In the following we 
shall consider only the range p G (1, oo). 

The link with the classical derivatives is given by this well-known fact: 
A function / belongs to W l ' p (Q) if and only if it is a.e. equal to a function 
u, absolutely continuous on almost all line segments in Q parallel to the 
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coordinate axes, whose (classical) derivatives belong to L p (f2) (the Beppo- 
Levi theorem). 

For 1 < p < oo and noninteger s > the Sobolev space W s,p (f2) of order 
s is defined as the space of all / with the finite norm 

( sr^ f f \D a f(x) - D a f(y)\ pX l/p 

WfWw^a) = [\\f\\ P w , ](n) + ^ JJ a \ x _ yl N +Pis - [s]) 

\a\ = [s] 

where [s] is the integer part of s. Details e.g. in [Adj . [Z] . 



2.2 Imbedding theorems 

One of the most useful and important features of the functions in Sobolev 
spaces is an improvement of their integrability properties and the compact- 
ness of various imbeddings. Theorems of this type were first proved by 
Sobolev and Kondrashev. Let us agree that the symbol > and ! -^ c — > stands 
for an imbedding and for a compact imbedding, respectively. 

Theorem 2.1. Let Q be a Lipschitz open. Then 

(i) If sp < N, then W S > P (Q) ^ L p * (f2) with p* = Np/(N - ps) (the 
Sobolev exponent). If\0\ < oo, then the target space is any L r (f2) with 
< r < p* . 

IfQ is bounded, then W S > P (Q) L q {Q) for all 1 < q < p* . 

(ii) If sp > N, then W j+S ' p (f2) C j (Q). If Q has the Lipschitz boundary, 
then W j+S > p (f2) ^ C^(77). 

If sp > N, then W j+S > p (f2) C j (Q) and W j+S ' P (Q) W*(n) 

for all 1 < q < oo. //, moreover, Q has the Lipschitz boundary, then 
the target space can be replaced by C^'^{Q) provided sp > N > (s — l)p 
and < /i < s — N/p. 

Note that if the imbedding W S,P {Q) ■=— > L q (Q) is compact for some q > p, 
then \ f2\ < oo. Moreover, if limsup r ,_ >00 \ {x G Q] r < \x\ < r + 1}| > 0, then 
W s,p (f2) L q (Q) cannot be compact. 
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2.3 Traces and Sobolev spaces of negative order 

Let s > and let Q be, for simplicity, a bounded open subset of M. N with 
boundary T of class C^'' 1 . Then with help of local coordinates, we can define 
Sobolev spaces W S,P (T) (also denoted H S (T) for p = 2) on T = df2 (see e.g. 
[N] , |Ad] for details). If / G C(f2), then /|r has sense. Introducing the space 
T>(f2) of restrictions in i? of functions in P(R ), one can show that if / G 
T>(f2), we have ||/|r||wi-Vp.p(r) — Cll/llw 1 ^^) so that, in view of the density 
of in iy 1,p (i7), the restriction of / to F can be uniquely extended to the 
whole W 1 ' P {Q). The result is the bounded trace operator 70 : W 1,p (f2) — > 
W /1_1//p,p (r). Moreover, every g G V(/ 1_1 /P'P(r) can be extended to a (non 
unique) function / G W /1,p (i7) and this extension operator is bounded with 
respect to the corresponding norms. 

More generally, let us suppose T is of class C* -1 ' 1 and define the operator 
Tr n for any / G V(Q) by Tr n / = (70/, 71/, . . . , 7fc-i/), where 

is the j-th order derivative of / with respect to the outer normal n at 
x G T; by density, this operator can be uniquely extended to a continu- 
ous linear mapping defined on the space W ,p (f2); moreover, ■jo(W k,p (f])) = 
W k ~ 1/p ' p (T). 

The kernel of this mapping is the space W k,p (f2) (denoted by Hq(Q) 
for p = 2), where W s,p (f2) is defined as the closure of T>(Q) in W S,P {Q) 
(s > 0). For 1 < p < 00 the following holds: W SJP (R N ) = W S ' P (R N ), 
W s >P(f2) = W s ' p ((2) provided < s < l/p. If s < 0, then the space W s > p (f2) 
is defined as the dual to W~ s,p (f2), where p' = p/(p — 1) (see e.g. |Trl] . 
|Tr2j ) . Observe that for an arbitrary Q a function / G W 1,p (fl) has the zero 
trace if and only if f{x)J dist(x, Y) belongs to L P (Q). 

For p = 2, we simply denote by H~ k (f2) the dual space of Hq(Q). In 
the case of bounded opens, we recall the following useful Poincare-Friedrichs 
inequality (for simplicity, we state it here in the Hilbert frame): 

Theorem 2.2. Let Q be bounded (at least in one direction of the space). 
Then there exists a positive constant Cp(f2) such that 

IMU a( fl) < C P {n)\\Vv\\ {L2mN for all v G H^H). (2.2) 
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2.4 The whole space case: Riesz and Bessel potentials 

The Riesz potentials X a naturally occur when one defines the formal pow- 
ers of the Laplace operator A. Namely, if / e S(M. N ) and a > 0, then 
T [(-A)"/ 2 /] (0 = l£| Q -77(0- This can be taken formally as a definition of 
the Riesz potential X a on S'(M. N ), 

i a f\.)=f- 1 [ier a ^/(o] (•) 

for any a Gi If < a < N, then I a f(x) = (I a * f)(x), where I a is the 
inverse Fourier transform of \C,\~ a , 

I a (x) = C a \x\ a ~ N , C a = Y ((N - a)/2) {i, N l 2 2 a Y{a/2))~ X 

(T is the Gamma function) is the Riesz kernel. The following formula is also 
true: 



IJx) = C n I t a 2 N e 



n\*\ 2 Olt 



o t 

Recall that every / e S(WL N ) can be represented as the Riesz potential I a g 
of a suitable function g E S(M. N ), namely, g = (— A) Q / 2 /; we get the repre- 
sentation formula 



f(x) =l a g(x) =C a 

Jk 



9iy) dy. 



N-a 



The standard density argument implies then an appropriate statement for 
functions in W k,p (M> N ) with an integer k and for the Bessel potential spaces 
H a,p (l Ar ) — see below for their definition. The original Sobolev imbedding 
theorem comes from the combination of this representation and the basic 
continuity property of I a , ap < N, 

1 1 rv 

I a : L?(R N ) - L q (R N ), - = - - ° 

q p JM 

To get an isomorphic representation of a Bessel potential space (of a 
Sobolev space with positive integer smoothness in particular) it is more con- 
venient to consider the Bessel potentials (of order a e M), 

Q a f(x) = (G a * /)(*) = ([1 + |e| 2 ]- Q/2 ^/(0) (x) 
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(with a slight abuse of the notations); the following formula for the Bessel 
kernel G a is well known: 

G a{x) = c a L I t 2 e t 4. 

Jo t 

(cf. the analogous formula for I a ), where c a = (Aii) a ^ 2 T(a/2). The kernels G a 
can alternatively be expressed with help of Bessel or Macdonald functions. 

Now we can define the Bessel potential spaces. For set and 1 < p < oo, 
let H S > P (R N ) be the space of all / e S'(R N ) with the finite norm 

\\f\\n^) = (7 f- 1 ((i+ \ey /2 rm)Y ^) 1/p . 

In other words, the spaces TC S,P (R ) are isomorphic copies of L P (R N ). 

For k = 0,1,2,..., plainly H k ' 2 (R N ) = W k ' 2 (R N ) by virtue of the 
Plancherel theorem. But it is true also for integer s and general 1 < p < oo 
(see e.g. |Trlj ). 

Remark 2.3. Much more comprehensive theory of general Besov and 
Lizorkin-Triebel spaces in M. N has been established in last decades, relying on 
the the Littlewood-Paley theory. Spaces on opens can be defined as restric- 
tions of functions in the corresponding space on the whole M. N , allowing to 
derive their properties from those valid for functions on R . The justification 
for that are extension theorems. In particular, there exists a universal exten- 
sion operator for the Lipschitz open, working for all the spaces mentioned up 
to now. We refer to [Trl] . JTr2] . 

2.5 Unbounded opens and weighted spaces 

The study of the elliptic problems in unbounded opens is usually carried out 
with use of suitable Sobolev weighted space. The Poisson equation 

- Ait = / in R N , N > 2, (2.3) 

is the typical example; the Poincare inequality ( 12.21) is not true here and it 
is suitable to introduce Sobolev spaces with weights. 

Let m E N, I < p < oo, a ER, k = m-N/p-a if N/p + a E {1, . . . ,m} 
and k — — 1 elsewhere. For an open Q C R N we define 

W™> P (H) = {v E V{Q), < |A| < k, p a -"H A l (log p)- 1 D x u E L P {Q), 
k + 1 < |A| < m, p a - m +\ x \D x u E L p ((2)}, 
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where p(x) = (1 + |x| 2 ) 1//2 , logp = log(2 + |x| 2 ). Note that W™ ,p is a reflexive 
Banach space for the norm ||.||w™- p defined by 

II*- = E \\P a - mMa \^ gP )- 1 D\\\l v{Q) + E \\p a - m+lal D X u\\ P LP ^y 

0<|A|<fc fc+l<|A|<m 

We also introduce the following seminorm 

i/p 



Mwc ,p 



(E iip Q ^ir^)) 



v |A|=m 

Let 

W™> p (ft) = {ve W™' p ; l0 {v) = . .. lm - x {v) = 0}. 

If ft is a Lipschitz domain, then W™> p (ft) is the closure of V(f2) in W™' p (ft), 
while V(ft) is dense in W™> p (ft). We denote by WI™' p \ft) the dual of 
W™' p (ft). We note that these spaces contain also polynomials, 



if f + a £ Z 




where [s] is the integer part of s and Pm = {0} if [s] < 0. The fundamen- 
tal property of functions belonging to these spaces is that they satisfy the 
Poincare weighted inequality. A open ft is an exterior domain if it is the 
complement of a closure of a bounded domain in R . 

Theorem 2.4. Suppose that ft is an exterior domain or ft = or ft = R . 
Then 

(i) the seminorm \ ■ \wa' p (o) ^ s a norms on W™' p (ft)/Pj, equivalent to the 
quotient norm with j' = min(m — 

(ii) the seminorm | • \w™' p (n) ^ s equivalent to the full norm on W™' p (ft). 



3 Variational approach 

Let us first describe the method on the model problem fll.ll) - fll.2p . supposing 
/ G L 2 (ft) and ft bounded. We first suppose that this problem admits a 
sufficiently smooth function u. Let v be any arbitrary (smooth) function; we 
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multiply equation (II. ip by v(x) and integrate with respect to x over i?; this 
gives 



-(Auv)(x) dx = / (fv)(x)dx. 
a J q 

Using the following Green's formula (dT(x) denotes the measure on T = dfl 

and = Vit( x) ■ n(x), where n(x) is the unit normal at point x of T 

oriented towards the exterior of i?) 

/ (Auv)(x)dx = - [ (Vu • Vv)(x) dx + [ {^v)(a)da, (3.1) 
Jn Jn Jr un 

we get, since v\r = 0: A{u,v) = L{y), where we have set 

A(u,v) = / Vm(x) • Vf (x) dx, L(v) = / f(x)v(x)dx. (3.2) 
J a J n 

The idea is to study in fact this new problem (showing first its equivalence 
with the boundary value problem), noting that it makes sense for far less 
regular functions u, v (and also /), in fact m,d6 Hl(Q) (and / G H~ l (Q)). 



3.1 The Lax-Milgram theorem 

The general form of a variational problem is 

to find u G V such that A{u, v) = L{y) for all v G V, (3.3) 

where V is a Hilbert space, A a bilinear continuous form defined on V x V 
and L a linear continuous form defined on V. We say moreover that A is 
V -elliptic if there exists a positive constant a such that 

A{u,u) > a\\u\\y for all u G V. (3.4) 

The following theorem is due to Lax and Milgram. 

Theorem 3.1. Let V be a Hilbert space. We suppose that A is a bilinear 
continuous form onVxV which is V -elliptic and that L is a linear continuous 
form on V . Then the variational problem ( 13.31) has a unique solution u on 
V . Moreover, if A is symmetric, u is characterized as the minimum value 
on V of the quadratic functional E defined by 

for all v G V, E{v) = -A{v, v) - L{v). (3.5) 
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Remark 3.2. (i) We have the following "energy estimate": \\u\\v < ll^llv- 
In the particular case of our model problem, this inequality shows the con- 
tinuity of the solution u G Hq(Q) with respect to the data / G L 2 (Q) (that 
can be weakened by choosing / G H^ 1 {Q)). 

(ii) Theorem 3.1 can be extended to sesquilinear continuous forms A defined 
on V x V; such form is called V^-elliptic if there exists a positive constant a 
such that 

ReA(u,u) > a\\u\\y for all u EV. (3.6) 
Finally, V is the dual of V. 

(iii) Denoting by A the linear operator defined on the space V by A(u, v) = 
(Au,v)v',v, for all v G V, the Lax-Milgram theorem shows that A is an iso- 
morphism from V onto its dual space V, and the problem (13.31) is equivalent 
to solving the equation Au = L. 

(iv) Let us do some remarks concerning the numerical aspects. First, this 
variational formulation is the starting point of the well known finite element 
method: The idea is to compute a solution of an approximate variational 
problem stated on a finite subspace of V (leading to the resolution of a 
linear system), with a precise control of the error with the exact solution u. 
Second, the equivalence with a minimization problem allows the use of other 
numerical algorithms. 

Let us now present some classical examples of second order elliptic prob- 
lems than can be solved with help of the variational theory. 

3.2 The Dirichlet problem for the Poisson equation 

We consider the problem on a bounded Lipschitz open Q C ~SL N 

—Au = f in i?, 

u = uq on 1 = ai/, 

with Uq G /f 1 / 2 (r), so that there exists Uq G H 1 ^) satisfying 7o(f^o) = u o- 
The variational formulation of problem (13.71) is 

to find u G Uq + Hl(ft) such that for all v G H^f}), A{u, v) = L(v), 

(3.8) 

with A given by (13.21) and a more general L with / G H x (f2), defined by 

L(v) = (f,v) H -i (a)tH i( a) . (3.9) 



11 



The existence and uniqueness of a solution of (13.81) follows from Theorem 3.6 
(and Poincare inequality (12. 2p . Conversely, thanks to the density of X>(i?) in 
Hq({1), we can show that u satisfies (13 .7p . More precisely, we get: 

Theorem 3.3. Let us suppose f E H~\Q) and u E H l / 2 {T); let U E 
H l (f2) satisfy 7q(Z7o) — u o- Then the boundary value problem (\3.7\i has a 
unique solution u such that u — UqE Hq(Q). This is also the unique solution 
of the variational problem ( 13.81) . Moreover, there exists a positive constant 
C = C{Q) such that 

\\u\\m{Q) < C + \\u \\ H i/2 {r) ) , (3.10) 

which shows that u depends continuously on the data f and uq. 

Moreover, using technics of Nirenberg's differential quotients, we have the 
following regularity result (see e.g. [Gri] ) : 

Theorem 3.4. Let us suppose that Q is a bounded open subset ofM, N with 
a boundary of class C 1 ' 1 and let satisfy f E L 2 (ft), uq E H 3 ^ 2 (T). Then 
u E H 2 (Q) and each equation in (13. 7p is satisfied almost everywhere (on Q 
for the first one and on T for the boundary condition) . Moreover, there exists 
a positive constant C = C(f2) such that 

IMIfl3(fl) < C [ll/IU 2 (fi) + |M|if3/2( r) ]. (3.11) 

By induction, if the data are more regular, i.e. / E H k (f2) and uq E 
#*+3/2(r) (with k E N), and if T is of class C k+1 '\ we get u E H k+2 (Q). 

Remark 3.5. Let us point out the importance of the open geometry. For 
example, if Q is a bounded plane polygon, one can find u E Hq{Q) with 
Am E C°°(i7), such that u H 1+1T l w (ft), where w is the biggest value of the 
interior angles of the polygon. In particular, if the polygon is not convex, 
the solution of the Dirichlet problem (13.71) cannot be in H 2 {Q). 

3.3 The Neumann problem for the Poisson equation 

We consider the problem (n is the unit outer normal on Y) 

—Am = / in i? 

du ^ (3.12) 

— — = h on 1 . 

on 
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Setting E(A) = {v G H l (Q); Av G L 2 (Q)}, the space V(Q) is a dense 
subspace, and we have the following Green formula for all u G E(A) and 
v G H l (Q): 

f f du 

I Au(x)v(x) dx = - / Vu(x) ■ Vv(x) dx + (—,70^)^-1/2^)^1/2^). 
Jq Jq On 

If u G if 1 (ft) satisfies (I3TT2D with / G L 2 (i?) and /i G if" 1 / 2 ^), then for any 
function v G H l (Q), we have, by virtue of the above Green formula, 



A(u, v) = L(v), Lv= (fv)(x) dx + (h, / yov) H -^(r),H 1 ^(r)- 



n 



But here, the form A is not if 1 (i7)-elliptic; in fact, one can check that, if 
problem f)3.12p has a solution, then we have necessarily (take v = 1 above) 

/ f(x)dx+ (h, l) H -i/2 (r)i ^i/2 (r) = 0. (3.13) 
J a 

Moreover, we note that if u is a solution, then u + C, where C is an arbitrary 
constant, is also a solution. So the variational problem is not well posed on 
H l (Q). It can be, however, solved in the quotient space if 1 (i7)/IR which is 
a Hilbert space for the quotient norm 

IH|tfi(fl)/R = inf ||u + fe||i,fl, (3.14) 

but also for the semi-norm v 1— > \v \m(n) — \/<A.(v, v), which is an equivalent 
norm on this quotient space, see [N] . 

Then, supposing that the data / and h satisfy the "compatibility condition" 
( 13.131) . we can apply the Lax-Milgram theorem to the variational problem 

to find u G H such that A(ii, v) = L(v) for all v G V (3.15) 

with V = if 1 (i7)/R. We get the following result (see e.g. [N]): 

Theorem 3.6. Let us suppose that Q is connected and that the data f G 
L 2 (i?) and h G H^^iY) satisfy (13. 13j) . Then the variational problem (13.151) 
has a unique solution u in the space if 1 (i7)/M and this solution is continuous 
with respect to the data, i.e. there exists a positive constant C = C(f2) such 
that 

\u\m{Q) < C (\\f\\i?(n) + \\h\\ H -i/2 {r) ) for all u G ii. 

Moreover, if V is of class C 1,1 and if the data satisfy / G L 2 (i7), g G 
if 1//2 (r), then every u G u is such that u G H 2 (Q) and it satisfies each 
equation in 13.121 almost everywhere. 
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3.4 Problem with mixed boundary conditions 

Here we consider more general boundary conditions: the Dirichlet conditions 
on a closed subset T\ of V — df2, and the Neumann, or more generally the 
"Robin", conditions on the other part r 2 = T — Ti. We seek u such that 
(/ G L 2 (f2), h G L 2 (r 2 ), a G L°°(r 2 )) 

—Am = / in f2, 

r + du h r (3 - 16) 

u = U on 1 1, au + — = h on 1 2 . 

on 

Let V = {v G 7 t> = on Tx}. Then ( 13. 3 p is the variational formu- 

lation of this problem with 

(i) A(u, v) = f n Vu(x) ■ Vv(x) dx + f r (a^y u 7 t> ) (a) da, 

(ii) L{y) = f n f(x)v(x) dx + f T {h^Qv){a) da. 

Supposing for example a > 0, we get a unique solution u G V for this 
variational problem by virtue of the Lax-Milgram theorem. Moreover, if 
u G H 2 (f2), then u is the unique solution in H 2 (Q) D V of the problem 

(M). 

3.5 The Newton problem for more general operators 

Let Q be a bounded open subset of K n . We now consider more general 
second order operators of the form v h- > — V.(MVv) + 6 • Vt> + cu, where 
b G [W 1 ' 00 ^)]* c G L°°(f2), M is an N x N square matrix with entries M ih 



and V • (MVf) stands for £) i7 - =1 /IM-^y J^l- We a ^ so assume that there i 



IS 



a positive constant % such that 

TV N 



i,j = l 1=1 



For given data / G L 2 (f2), h G £ 2 (r), we look for u solution of the problem 
-V ■ (MVu) + b ■ Vu + cu = f in Q, 

(3-17) 

au + n ■ (MVm) = h on T. 
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We assume that a G L°°(T). The variational formulation of this problem is 
still with V = H\Q) and 

A(u,v) = I MVu -Vv dx + / [b ■ Vu + cu]v dx + / a^ou^ov da, (3.18) 
if? Jn Jt 

L(v)= / f(x)v(x) dx + f(hj v)(a) da, (3.19) 
Jn Jt 

If the following conditions 

1 1 
c — -V ■ b > Co > a.e. on i?, a + -b • vgeC\ > a.e. on T 

are fulfilled, with (C , Ci) 7^ (0, 0), then the bilinear form A is V-elliptic and 
the Lax-Milgram theorem applies. 

3.6 A biharmonic problem 

We consider the Dirichlet problem for the operator of fourth order: (c G 
L°°(Q)): 

A 2 u + cu = f in 0, (3.20) 

u = u on T, — — — h on T. (3.21) 
on 

Theorem 3.7. Let us suppose that i? has a boundary of class C 1,1 and that 
the data satisfy f G H~ 2 {f2), u G H 3 / 2 (T), h G H l ' 2 (T). Let U G H 2 (Q) 
be such that 7o(?7o) = Uq,1x{Uq) = h. Then, if c > a.e. in Q, the boundary 
value problem f l3.20l) - fl3.21l) has a unique solution u such thatu—Uo G H^f}), 
and u is also the unique solution of the variational problem 

to find ueU + H 2 (f2) such that A{u, v) = l{v) for all v G H 2 ((2), (3.22) 

where l(v) = (f,v} H -i {n)>H 2 {n) and 

A(u,v) = / Au(x)Av(x) dx + / (cuv)(x) dx. (3.23) 
Jn Jn 

Moreover, there exists a positive constant C = C(f2) such that 

\\u\\ H 2 (n) < C[\\f\\ H -2 {n) + ||uo|| H 3/2 (r) + || h || #1/2 ( r)], (3.24) 
which shows that u depends continuously upon the data f , Uq and h. 
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Remark 3.8. The Hilbert space choice V is of crucial importance for the 
l^-ellipticity. In fact, let us consider for example the problem (I3.20p . (I3.25p . 
with 

dAu 

Au = on T, — — = on T. 3.25 
on 

In fact, the associated bilinear form is not ^-elliptic for V = H 2 (f2) but it 
is V-effiptic for V = {v G L 2 (f2); Av G L 2 Q)}. 



4 General elliptic problems 

Here £2 will be a bounded and sufficiently regular open subset of R N . Let us 
consider a general linear differential operator of the form 

A(x, D)u = a^D^u, a^x) G C (4.1) 
\tA<i 

Setting A (x, £) = Y^\n\=i a /j.( x )C tl j we sa y that the operator A is elliptic at 
a point x if A (x, f) ^ for all £ eR N - {0}. One can show that, if iV > 3, 
/ is even, i.e. / = 2m; the same result holds for iV = 2 if the coefficients a M 
are real. Moreover, for N > 3, every elliptic operator is properly elliptic, in 
the following sense: For any independent vectors £, £' in M. N , the polynomial 
t i — > Aq(.,C, + r£') has m roots with positive imaginary part. 

The aim here is to study boundary value problems of the following type: 

Au = f in f], (4.2) 

BjU = gj on T, j — 0, ...,m — 1, (4.3) 

where A is properly elliptic on i?, with sufficiently regular coefficients, and 
the operators Bj are boundary operators, of order mj < 2m — 1, that must 
satisfy some compatibility conditions with respect to the operator A (see 
|RR| for details; these conditions were introduced by Agmon, Douglas and 
Nirenberg). For example, A = (— l) m A m and Bj = j^j is a convenient choice. 

In order to show that problem (14. 2 1) - 04. 3 p has a solution u G H 2m+r {f}) 
(r G N), the idea is to show that the operator V defined by u i— ► V(u) = 
(Au, Bqu, B m ^iu) is an index operator from H 2m+r (f2) into G = H r (f2) x 
Ti™~^H 2m+r ~ m: i~ l l 2 (T) and to express the compatibility conditions through 
the adjoint problem. 

We recall that a linear continuous operator V is an index operator if 
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(i) dim Ker V < oo, and Im V closed 

(ii) codim Im V < oo. 

Then the index is given by x(^) = dim Ker P - codim Im V. We 
recall the following Peetre's theorem: 

Theorem 4.1. Let E, F and G be three reflexive Banach spaces such that 
E F , and V a linear continuous operator from E to G. Then condition 
(i) is equivalent to 

(Hi) there exists C > 0, such that for all u G E, we have 

\\u\\e < C {\\Vu\\ G +\\u\\ F .) 

Applying this theorem to our problem (I4.2p - (I4.3I) . condition (i) results 
from a priori estimates of the following type: 

\\u\\ H 2m + r^ < C (||"Pm|| G + \\u\\ H 2m + r-l^) 

and condition (ii) by similar a priori estimates for the dual problem. 

5 Second order elliptic problems 

We consider a second order differential operator of the "divergence form" 

N N 

Au = - (x)u Xi ) Xj + Y b \ x ) u ^ + c{x)u (5.1) 

i,j=X i=l 

with given coefficient functions a w , b l , c (i, j = 1, . . . , N), and where we have 
used the notation u Xi = J^-. Such operators are said uniformly strongly 
elliptic in i? if there exists a > such that 

"' J (-'-K'C > a|£| 2 for all x 6 f2, ^ G R N . 

|i|=lil=l 

Remark 5.1. There exists elliptic problems for which the associated vari- 
ational problem does not necessarily satisfy the ellipticity condition. Let 
us consider the following example, due to Seeley: Let Q = {(r, 6) G 
(rr,2n) x [0,2tt]} and A = -(e ie ^) 2 - e 2i %l + fy. One can check that, 
for all A G C, the problem Au + Xu = f in i? and u = on T admits nonzero 
solutions u which are given by (with [i such that fi 2 = A) u = sinr cos(fie~ t6 ) 
and u = sin rsin(fie~ te ) for A ^ 0; « = sinr and u = sm8e~ ie for A = 0. 
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Most of results concerning existence, unicity, regularity for second order 
elliptic problems can be established thanks to a maximum principle. There 
exist different types of maximum principles, that we now present. 

5.1 Maximum principle 

Theorem 5.2 (weak maximum principle). Let A be a uniformly strongly 
elliptic operator of the form ( 15. ip in a bounded open Q C R , with a % \ b l , c G 
L°°(f2) and c > 0. Let u G C 2 (^) n C(?2) and 

An > [resp. Au < 0] in fi. 

Then 

infw > infii - [resp. supw < supit + ], 

where u + = max(w, 0) and u~ = — min(w,0). If c = in Q, one can replace 
u~ [resp. u + ] by u. 

Theorem 5.3 (strong principle maximum). Under the assumptions of the 
above theorem, if u is not a constant function in C 2 (Q) fi C(i?) such that 
Au > [ resp. Au < 0] , then inf^ti < u(x) [resp. sup^w > u{x)], for all 
x G Q. 

Remark 5.4. These two maximum principles can be adapted to elliptic 
operators in nondivergence form, i.e., 

N N 

AU = - " '(■>') 11 -r,.r, + Y b \ X ) U *i + C ( X ) U - ( 5 - 2 ) 

i,j=l i=l 

5.2 Predholm alternative 

We now present some existence results which are based on on the Fredholm 
alternative rather than on the variational method. 

Let us consider two Hilbert spaces V and H, where V is a dense subspace 
of H and V H . Denoting by V the dual space of V, and identifying H 
with its dual space, we have the following imbeddings: V > H °-> V. Let 
A be a sesquilinear form on V x V, ^-coercive with respect to H, that is, 
there exists A G R and a > such that 

Re(A(v, v)) + Xo\\ v \\h — a \\ v \\v ^ or a ^ v eV . 
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Denoting by A the operator associated with the bilinear form A (see item 
(hi) of Remark 3.1), the equation Au = f is equivalent to u — \ Tu = g, 
with T = (A + Xold)' 1 and g = Tf. Note that T is an isomorphism from H 
onto D{A) = {u E H; Au E H}). 

The operator T : H — > if is compact and, thanks to the Fredholm alter- 
native, there are two situations: 

(i) either Ker A = and A is an isomorphism from D(A) onto H 

(ii) or Ker A^O; then Ker A is of finite dimension, and the problem Au = f 
with / G H admits a solution if and only if / G lm.A = [Ker(A*)] ± . 

We now give another example in a non Hilbertian frame. Let us consider 
the problem |Grij : Au = f in Q and Bu = g on T where T is of class C 1 ' 1 , 
A, which is defined by (15.11) . is uniformly strongly elliptic with a*- 7 = a? % G 
C 0,1 (J?), 6\c G L°°(i?) and £?u = 7o(«) or Bu = ji(u). One can show 
that the operator u t— > {Au, Bu) is a Fredholm operator of index zero from 
W 2 > p (f2) in L p {Q) x jy^-Vp.pjT) (with d = if Bu = 7 (u) and d = 1 if 
Bu = 7i (it)). 



5.3 Regularity 

Assume that i? is a bounded open. Suppose that u G Hq{Q) is a weak 
solution of the equation 

Au = f in i?, 

it = on 1 , 

where A has the divergence form (15. ip . We now address the question whether 
u is in fact smooth: this is the regularity problem for weak solutions. 

Theorem 5.5 (if 2 -regularity) . Let Q be open, of class C 1 ' 1 , a % i G C l (Q), 
b l ,c G L°°(f2), f G L 2 (Q). Suppose furthermore that u G ii 1 (i7) is a weak 
solution of (15. 3p . Then u G if 2 (i?) and we have the estimate 

\\u\\ H *{n) < C(\\f\\ L 2 { n) + \\u\\ L 2 {n) ), 

where the constant C depends only on Q and on the coefficients of A. 

Theorem 5.6 (higher regularity). Let m be a nonnegative integer, Q be 
open, of class C m+1 ^ and assume that a ij G C m+l (j2), b\c G C m+1 (Q), 
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/ G H m (f2). Suppose furthermore that u G H l (Q) is a weak solution of 
(jCT) . T/ien it G H m+2 (Q) and 

||w||i?™+ 2 (fi) < c(||/||jf m (fi) + II^IIl 2 (^)), 

where the constant C depends only on Q and on the coefficients of A. In 
particular, if m > N/2, then u G C 2 (f2). Moreover if Q is of C°° class and 
f G C°°(77) ; a'i G C°°(n), b\ce C°°(72), then u G C°°(72). 

Remark 5.7. (i) If u G Hq({2) is the unique solution of (15.31) . one can omit 
the L 2 norm of u in the right hand side of the above estimate, 
(ii) Moreover, let us suppose the coefficients a t; i , h l and c are all C°° and 
/ G C°°(i7); then, if u G H l {Q) satisfies = /, u G this is due to 

the "hypoellipticity" property satisfied by the operator A. 

We have a similar result in the LP frame [Gri] : 

Theorem 5.8 (W /2 ' p -regularity). Let Q be open, of class C 1,1 , a lj G C 1 (J7) ; 
b l ,c G L°°(Q). Suppose furthermore that b % = 0, 1 < i < N and c > a.e 
or c > (3 > a.e . Then for every f G L p (f2) there exists a unique solution 
u G W 2 ' p (f2) of (ED- 

6 Unbounded open 
6.1 The whole space 

Note in passing that we shall work with the weighted Sobolev spaces W™ ,p (f2) 
defined in Subsection 12.51 

Theorem 6.1. The following claims hold true: 

(i) Let f G Wq 1,p (R n ) satisfy the compatibility condition 

(f, 1 ) w - 1 ^( R iV) xWo i.p' (R iv ) =0 if p >N. 

Then the problem (12. 3p has a solution u G W Q ' P (M. N ), which is unique 
up to an element in "P[i_jv/pi an d satisfies the estimate 

\\ U \\w^ p (R N )/V {1 ^ N/p] - CWfWw- 1 '"^)- 

Moreover, if 1 < p < N, then u = E * f . 
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(ii) If f e L P (R N ), then the problem (J23D has a solution u G Wq' p (R n ), 
which is unique up to an element in V\2-N/p] an d ifl<p< N/2, then 
u — E * f. 

The Calderon-Zygmund inequality 

d 2 <p 



LP 



and Theorem 12.41 are crucial for establishing Theorem 16 . 11 Further, point (i) 
means that the Riesz potential of second order satisfies 

h : W^' p (R N )±V^ N/pf] - W^{R N )/V^ 



(where the initial space is the orthogonal complement of V[i-n/p'] in 
W^ P {R N )) and it is an isomorphism. 
Note that here 

Wo' p (R N ) = {v G L P *(R N ); Vv G L P (R N )} 

for 1 < p < N and 1/p* = 1/p — 1/N. And for 1 < r < y, we have also the 
continuity property 

I 2 : L r (R N ) -> L q (R N ), for - = 

q r N 

Remark 6.2. The problem 

u-Au = f mR N (6.1) 

is of a completely different nature than the problem (12. 3p . The class of 
function spaces appropriate for the problem (16. ip are the classical Sobolev 
spaces. With help of the Calderon-Zygmund theory one can prove that if / G 
L P (R N ), then the unique solution of (16.11) belongs to W 2 ' P (R N ) and can be 
represented as the Bessel potential of second order (see [S]): u = G* f, where 
G is the appropriate Bessel kernel, that is, G, for which G(£) ~ (1 + |^| 2 ) -1 ^ 2 - 
Recall that in particular G{x) ~ |a;|~ 1 e - 'l x ' for N = ?>. In the Hilbert case 
/ G L 2 (R N ) we get 

(l + \C\ 2 )ueL 2 (R N ), 

which, by Plancherel's theorem, implies that u G H 2 (R N ). For / G 
W~ 1,P (R N ), the problem 16.11 has a unique solution u G W l,p {R N ) satisfy- 
ing the estimate 

IMIwi.p(RJV) < C(p,n) || /|| w -i, P r 
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6.2 Exterior domain 



We consider the problem in an exterior domain with the Dirichlet boundary 
condition 

—Au = f in 42, . 
u = g on OiJ, 

where / G W / " 1 ' p (i2) and 3 G ^"^(di?). Invoking the results for R N 
and bounded domains one can prove the existence of a solution u G W ' p (f2) 
which is unique up to an element of the kernel Aq(Q) = {z6 Wq' p {Q); Az = 
0} provided that / satisfies the compatibility condition 

(M = (g,^) for all (p G Aq (f2). 

The kernel can be characterized in the following way: It is reduced to {0} if 
p = 2 or p < N and if not, then 

A p (f2) = {C{\ - 1); C G R} if p > N > 3, 

where A is (unique) solution in W ' (i?) fl Wo' p (i2) of the problem AA = 
in i2 and A = 1 on df2, and 

Ag(X2) = - M ); C G M} if p > iV = 2, 

where «o(^) = (27r | T | ) — 1 J r log \y — x\da y and /i is the only solution in 
Wq' 2 ((2) n W 1,P (J?) of the problem A/i = in Q and fi = u on T. 

Remark 6.3. Similar results exist for the Neumann problem in an exterior 
domain (see |AGGJ). The framework of the spaces VK™' P (R+) also for the 
Dirichlet problem in was considered in the literature, too. For a more 
general theory see [KM . 



7 Elliptic systems 
7.1 The Stokes system 

The Stokes problem is a classical example in the fluid mechanics. This system 
models the slow motion with the field of the velocity u and the pressure 7T, 
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satisfying 

— vAu + 'Vir = f in f2, 

divu — h in Q, (S) 

u = g on T = df2, 

where v > denotes the viscosity, / is an exterior force, g is the velocity of 
the fluid on the domain boundary and h measures the compressibility of the 
fluids (if h — 0, it is an incompressible fluid). The functions h and g must 
satisfy the compatibility condition 

h(x) dx = J g-nda. (7-1) 

Theorem 7.1. Let Q be a Lipschitz bounded domain in ~R N , N > 2. Let 
f G H-\G) N , h G L 2 (J?) and g G H l ' 2 {V) N satisfy (1711) . T/ien tte pro6/em 
(S*) /ias a unique solution (u, n) G if 1 (i7) 7V x L 2 (i7)/IR satisfying the a priori 
estimate 

INI-ffi(fl) + IKIIl 2 (/2)/k < cdl/llff- 1 ^) + II^IIl 2 («) + 11^ II jyi/2(r) ) • 

In order to prove Theorem 17.11 one can start with a homogeneous problem. 
The procedure of finding u is a simple application of the Lax-Milgram theo- 
rem. Application of De Rham's theorem gives the pressure tt. We introduce 
the space 

V = {ve V{Q) N - div{7= 0} 
and define F G H-\n) N by 

{F,v) H -i X ffi = for all v G V. 

Moreover, there exists tt G L 2 (f2), unique up to an additive constant, and 
such that F = Vtt. The problem (S), which we transform to the homoge- 
neous case (h = 0, g = 0), can be formulated on an abstract level. Let X and 
M be two real Hilbert spaces and consider the following variational problem: 
Given L G X' and X G M', find (u,tt) <E X x M such that 

A(u,v) +B[v,ir) = L(v), v e X, 
B[u,q]=X(q), qeM. 

where the bilinear forms A, B and the linear form L are defined by 

A(u,v) = / W-W, 5[i?,g] = - / [qV ■ v\, L(v) = f ■ v. 
Jn J q Jn 
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Theorem 7.2. If the bilinear form A is coercive in the space 

V = {vEX; B[v, q] = for all q G M, 

i.e. if there exists a > such that 

A(v,v) > a\\v\\x, v E V, 

then the problem (17.21) has a unique solution (u, it) if and only if the bilinear 
form B satisfies the "inf-sup" condition: 

B(v, q) 

there exists (3 > such that inf sup -rr-^r — ^—r — > (3. 

qeM# eX \\v\\x\\q\\M 

As for the Dirichlet problem, the regularity result is the following: 

Theorem 7.3. Let Q be a bounded domain in M. N , of the class C m+1 ' 1 if 
m G N and C 1 ' 1 if m = -1. Let f G W m > p {Q) N , h G W m+1 ' p ((2) and 
g G W m+2 ~ 1 / p ' p (T) N satisfy condition ( 17. ip . Then the problem (S) has a 
unique solution (u,ir) G W m+2 ' P ((2) N x W m+1 > p (f2)/R. 

Remark 7.4. It is possible to solve (S) under weaker assumption, for in- 
stance if f e W-^ p (n'), h = and g G ^^/"(r)^. We can prove that 
then (u,ir) G L p (Q) n x W~ l ' p (Q). 

7.2 The linearized elasticity 

The equations governing the displacement u— (ui, 112,113) of a three dimen- 
sional structure subjected to an external force field / are written as (Q is a 
bounded open subset of M 3 and Y = dft) 

-fiAu - (A + /i)V(V • u) = f mf2, 

u — on r , 

3 

^2 = (ji on T x = r - r , 

3=1 

where A > and /i > are two material characteristic constants, called the 
Lame coefficients, and (v = (t>i, t> 2 , t> 3 )) 

3 

aij(v) = aji(v) = XSy ^e kk {v) + 2y£ ij {v) 

k=i (7.3) 

with £ij(v) = Eji(v) = -{djVi + diVj), 
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where Sij denotes the Kronecker symbol, i.e. 5^ = 1, for i = j and 5ij = 0, for 
i ^ j. These equations describe the equilibrium of an elastic homogeneous 
isotropic body that cannot move along r ; along Ti, surface forces of density 
g = (gi, #2j gs) are given. The case T± = physically corresponds to clamped 
structures. The matrix with entries £ij{u) is the linearized strain tensor while 
<Jij{u) represents the linearized stress tensor; the relationship (17.31) between 
these tensors is known as Hooke's law. We refer for example to [CLj . jNH] 
(and to the references herein) for most of the results stated in this paragraph. 
The variational formulation of this problem is 

to find u G V such that A(u, v) = L(v) for all v £ V, (7-4) 

where the bilinear form A and the linear form L are given by 




L(v) = / f(x) ■ v(x) dx + g(x) ■ v{a) da; (7.5b) 

J Q JFx 

The functional space V is defined as 



V = {v = (v u v 2 ,v 3 ) e [H x {a)f- lQVi = onr , 1 < i < 3}. 



To prove the ellipticity of A, one needs the following Korn inequality: There 
exists a positive constant C(i?) such that, for all v = (vi,v 2 ,v 3 ) G [if 1 (i?)] 3 , 
we have 



Hli,«<c(fi) 



E 

.«,j=i 



\ £ ij{ V )\\L2(n) 



\L 2 {Q) 



i=l 



1/2 



(7.6) 



The following result holds true: 

Theorem 7.5. Let Q be a bounded open in IR 3 with a Lipschitz boundary, 
and let To be a measurable subset of V , whose measure (with respect to the 
surface measure dT{x)) is positive. Then the mapping 



■nl/2 



E 



\ £ iA v )\\L^(Q) 



is a norm on V , equivalent to the usual norm || . 
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As a consequence, we get: 

Theorem 7.6. Under the above assumptions, there exists a unique u G V 
solving the variational problem (17.41) - 07. 5p . This solution is also the unique 
one which minimizes the energy functional 

3 

m = \ [ [\{V-vf + 2^[e l3 {v)f]{x)dx 
— / fix) ■ v(x) dx + / g(x) ■ v(a) da 

over the space V . 
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